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Abstract 

An interesting discovery in the last two years in the field of mathematical physics has been the 
exceptional Xg Laguerre and Jacobi polynomials. Unlike the well-known classical orthogonal poly- 
nomials which start with constant terms, these new polynomials have lowest degree t = 1,2,..., 
and yet they form complete set with respect to some positive-definite measure. While the mathe- 
matical properties of these new Xg polynomials deserve further analysis, it is also of interest to see 
if they play any role in physical systems. In this paper we indicate some physical models in which 
these new polynomials appear as the main part of the eigenfunctions. The systems we consider 
include the Dirac equations coupled minimally and non-minimally with some external fields, and 
the Fokker-Planck equations. The systems presented here have enlarged the number of exactly 
solvable physical systems known so far. 



1 



I. INTRODUCTION 



The past two years have witnessed some interesting developments in the area of exactly 
solvable models in quantum mechanics: the number of exactly solvable shape-invariant 
models has been greatly increased owing to the discovery of new types of orthogonal poly- 
nomials, called the exceptional Xi polynomials. Two families of such polynomials, namely, 
the Laguerre- and Jacobi-type X\ polynomials, corresponding to £ = 1, were first proposed 
by Gomez-Ullate et al. in within the Sturm-Liouville theory, as solutions of second-order 
eigenvalue equations with rational coefficients. Unlike the classical orthogonal polynomials, 
these new polynomials have the remarkable properties that they still form complete set 
with respect to some positive-definite measure, although they start with a linear polyno- 
mials instead of a constant. The results in jl| were then reformulated by Quesne in the 
framework of quantum mechanics and shape-invariant potentials, first in 2J by the point 
canonical transformation method, and then in jsj by super symmetric (SUSY) method 4] 
(or the Darboux-Crum transformation 1 51). Soon after these works, such kind of exceptional 

n 

polynomials were generalized by Odake and Sasaki to all integral £ = 1,2,... [6| (the case 
of £ = 2 was also discussed in [3|). By construction these new polynomials satisfy the 
Schrodinger equation and yet they start at degree £ > instead of the degree zero constant 
term. Thus they are not constrained by Bochner's theorem {jj], which states that the orthog- 
onal polynomials (starting with degree 0) satisfying a second order differential equations can 
only be the classical orthogonal polynomials, i.e., the Hermite, Laguerre, Jacobi and Bessel 
polynomials. 

Later, equivalent but much simpler looking forms of the Laguerre- and Jacobi-type Xi 



polynomials than those originally presented in |6| were given in [8]. These nice forms were 
derived based on an analysis of the second order differential equations for the X^ polynomials 
within the framework of the Fuchsian differential equations in the entire complex x-plane. 
They allow us to study in-depth some important properties of the Xt polynomials, such 
as the actions of the forward and backward shift operators on the Xg polynomials, Gram- 
Schmidt orthonormalization for the algebraic construction of the X^ polynomials, Rodrigues 
formulas, and the generating functions of these new polynomials. 
Recently in 0] the X\ Laguerre polynomials in 



1[ are generalized to the Xi 
polynomials with higher £ based on the Darboux-Crum transformation. Then in 



iaguerre 



1 Of ] such 
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transformation was successfully employed to generate the Xe Jabobi as well as the Xe La- 
guerre polynomials as given in jg]. 

While the mathematical properties of these new polynomials deserve further analysis, 
it is also of interest to see if they play any role in physical systems. It is the purpose of the 
present work to indicate some physical models that involve these new polynomials. We shall 
be mainly concerned with the Xe Laguerre polynomials for clarity of presentation. Models 
that are linked with the X? Jacobi polynomials will be briefly mentioned at the end. 

The plan of the paper is as follows. First we review the deformed radial oscillators 
associated with the exceptional Xe Laguerre polynomials in Sect. II. In Sect. Ill the Dirac 
equation minimally coupled with an external magnetic field is considered. We present the 
forms of the vector potentials such that the eigenfunctions of the Dirac equation are related 
to the Xe Laguerre polynomials. Dirac equations with non-minimal couplings are then 
mentioned in Sect. IV, and the Fokker-Planck equations are considered in Sect. V. Sect. VI 
concludes the paper. 



II. EXCEPTIONAL X e LAGUERRE POLYNOMIALS 

Consider a generic one- dimensional quantum mechanical system described by a Hamilto- 
nian H = —d 2 /dx 2 + Vq(x). Suppose the ground state is given by the wave function 4>o(x) 
with zero energy: H(po = 0. By the well known oscillation theorem </>o is nodeless, and thus 
can be written as 0o = e w °^ x \ where Wq(x) is a regular function of x. This implies that the 
function W (x) completely determines the potential V : V = W$ + W$ (the prime here 
denotes derivatives with respect to x). Thus Wq{x) is sometimes called a prepotential. Con- 
sequently, the Hamiltonian can be factorized as H = Hq + ^ = AqAq, with = ±d/dx — WQ. 
It is trivial to verify Aq<J)q(x) = 0. 

It is a remarkable fact that most well-known one- dimensional exactly solvable systems 



possess a property called shape invariance ll|. This means the Hamiltonian Hq~\ de- 
fined by Hq~^ = AqAq with potential Wq — Wq, is related to Hq^ by the relation 
Hq \X) = Hq + \\ + 6) + £i(A). Here A = (Ai, A 2 , . . .) is a set of parameters of the Hamil- 
tonian H = H(X), 8 is a certain shift of these parameters, and £i(A) is some function of 
A. The mapping that relates the potentials and Hq~^ is called the Darboux-Crum 

transformation |5( {H^ and Hq ' are also called SUSY partners in the context of SUSY 
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quantum mechanics [4]). Shape invariance is a sufficient condition [12] that enables one to 
determine the eigenvalues and the corresponding eigenfunctions of Hq + \\) exactly (see 4] 
for details). Specifically, we have 

n-l 

£o(A) = 0, ^ m (A) = ^f 1 (A + A; < 5), (1) 
m (x;A) oc A-{\)A {\ + 6)---A~{\+{m-l)d)xe Woix ' x+mS \ m = l,2,.... (2) 

The new shape-invariant systems related to the exceptional polynomials are determined 
by certain prepotentials We (I = 1,2,...), which are obtained by deforming some shape 
invariant prepotentials Wq The i = case corresponds to the original system. Three 
families of such exactly solvable deformed systems were presented in |6j, which correspond 
to deforming the radial oscillator and the trigonometric/hyperbolic Darboux-Poschl- Teller 
potentials in terms of their respective eigenfunctions. 

For the purpose of this paper, we shall only consider the exceptional Xe Laguerre polyno- 
mials, which appear in the deformed radial oscillator potentials. The original radial oscillator 
potential is generated by the prepotential 

2 

COX 

W (x;g) = — + g\ogx, < x < oo. (3) 

Here A = g > and co > 0. The Hamiltonian is 



H { + \g) = A-(g)A+(g) = ~ + u 2 x 2 + - (2g + l)u>. (4) 



This potential is shape-invariant with shift parameter 8 = 1 and £\(g) = 4w: \g) = 
H^ + \g + 1) + 4u>. The eigen-energies and eigenfunctions are (n = 0, 1, 2, . . .) 

£ n {g) = Ancu; n (x; g) = e Wo ^P n { m g), P n ( V ; g) = L^' h) (v), (5) 

where rj(x) = uox 2 is one of the so-called sinusoidal coordinates fl3 |. 

Below we shall present radial oscillators related to the Laguerre polynomials. We first 
treat them as systems deformed from the original radial oscillator by appropriate deforming 
functions. Then we show how they can also be considered as the Darboux-Crum (or SUSY) 
partner of the original radial oscillator. 
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A. Deformed radial oscillators 



The Hamiltonian H^ + \g) of the deformed radial oscillator is 



H ( t\g) = ~^2+ W?(x;g) + W?(x;g), (6) 



where Wi is given by 



We(x; g) = -^ + (g + £) logx + log ^'f + 1 \ (7) 

Here £,e(r); g) is a deforming function. It turns out there are two possible sets of deform- 
ing functions ^i(r);g), thus giving rise to two sets of infinitely many exceptional Laguerre 
polynomials, termed LI and L2 type j^, 8]. These & are given by 

f L^i-n) :L1 
[ 4 9 %) : L2. 

For both types of the eigen-energies are £^ n (<?) = £n(g + £) = 4na;, which are indepen- 
dent of g and £. Hence the deformed radial oscillator is iso-spectral to the ordinary radial 
oscillator. The eigenfunctions are given by 

e ~\bJX 2 x g+l 

4>e,n(x; g) = — — r-Pi, n (ri\ g), (9) 

&(w g) 

where the corresponding exceptional Laguerre polynomials Pi,n(jl\9) {£ — 1)2,..., n = 
0, 1,2,.. .) can be expressed as a bilinear form of the original Laguerre polynomials and the 
deforming polynomials, as given in [8J]: 



Urn g + i)P n {w g + £-i)- Uv; g)d v P n {m g + £-i) = li 
PiAv;g)= < (n + g + \)- l {{g + \)amg + ^Pn(mg + z + ^) ( 10 ) 
[ +vUv;g)d v Pn(ri;g+t + i)) =L2. 

The Xg polynomials Pe, n (r); g) are degree £ + n polynomials in 77 and start at degree £: 
Pe,o(V'i g) = fLiiW-i g + !)■ They are orthogonal with respect to certain weight functions, which 
are deformations of the weight function for the Laguerre polynomials (for details, see js]). 



B. Darboux-Crum pairs 

As mentioned in Sect. 1, the new exceptional orthogonal polynomials have recently been 
re-derived from the corresponding ordinary polynomials based on the Darboux-Crum trans- 
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formation (9), [loj]. Below we shall reconsider the case of the deformed radial oscillator, 
following the approach of [lo| . 

Instead of finding the prepotential We that gives Hq + ^ as the deformed oscillator, as was 



done in the previous subsection, one determines a new prepotential Wi so that it is Hq ' 



that gives the deformed oscillator whi 
The new prepotential Wi(x\ g) are 



e Hq is related to the ordinary oscillator. 
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\ux 2 + (g + l-l)\ogx + \ogUv(x);g)i ^ > 1/2, : LI 
Wt(x;g) = { (11) 

■\ux 2 - {g + £)logx + logCMx);g), g > -1/2, : L2 

where ^e(rj;g) are as given in (jSJ). Now consider the pairs of Hamiltonians %g (#) and 
ft^fo) (£=1,2,...) defined by 

= ^(5)^(3), H { -\g) = Al(g)Ai(g), (12) 
^0?) = ±-^-Wi{x;g). (13) 

Using the differential equation for the Laguerre polynomial, the Hamiltonians H\ + \g) can 
be shown to be related to that of the radial oscillator HQ + \g) in as 

w , + , fe)H ^ + '( 9 + *-l) + 2(2 S + «-l). :L1 
HnS + f + l)+2(2 S + l) W :L2. 

The partner Hamiltonians are found to equal to the Sasaki-Odake Hamiltonian H^ + \g) in 
dSJ) up to additive constants: 

H { l + \g) + 2{2g + M- l)u : LI 
i/j +) (^) + 2(2^ + l)w :L2. 

10| that the partner Hamiltonian T-L^~\g) are exactly iso-spectral to the 



«r , fa) = <T, u '™ T ~ T (15) 



It is shown in 



radial oscillator Hamiltonian "H^ (<?), which have the following eigenvalues and the corre- 
sponding eigenfunctions (r] = ux 2 ): 

LI: < + n ) (9)=4(« + ff + 2f-^ Wj 0«(x^) = e-^V^- 1 Li 9+ '- f) (r ? ), (16) 

L2: = 4 L + g + ^ u, <f>$(x; g) = e^V^ 1 !^^). (17) 

The eigenfunctions (f)[~J(x;g) of l-if\g) are obtained by applying the operator Af(g) on 

^^(^S 9)'- 4>tn(. x 'i9) = ^{g)(jy^l{ x \ 9)- ^ turns out tnat 4>tn( x 'i9)-> U P to multiplicative 
constants, are just the eigenfunctions 4>e t n(vi9) m & an d ( ITOj) of the Sasaki-Odake Hamil- 
tonian. 
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III. DIRAC EQUATION WITH MAGNETIC FIELD IN CYLINDRICAL COOR- 
DINATES 



Having reviewed the deformed radial oscillators associated with the exceptional X e La- 
guerre polynomials, we will like to see if there are physical systems in which these new 
polynomials could play a role. It turns out that there are indeed such systems. In what 
follows we shall indicate some of them. 

First let us consider the Dirac equation in 2+1 dimensions coupling minimally with a 
cylindrically symmetric magnetic field. The discussion can be extended to 3+1 dimensions 
where the magnetic field does not depend on the variable z. The Dirac equation for a unit 
positive (q — 1) charged particle minimally coupled to an magnetic field B = V x A has 
the form (we set c = h = 1) 

H d V(t) = EV(r) 
H D = a • (p - A) + a 3 M, (18) 

where Hp is the Dirac Hamiltonian, p = — iV is the momentum operator, E and M are 
energy and rest mass of the particle, and a are the Pauli matrices. 

We shall consider magnetic field which is cylindrically symmetric. Then the vector po- 
tential has only 0-component: 

A(r) = A <t ,(r)<j>, r= \r\. (19) 

Instead of the variable x, we shall use the conventional notation r for the radial variable 
here, and in Sect. IV. A and B. The magnetic field is 

Bs(r) = ~(rMr))- (20) 
The wave function is taken to have the form 

V'm r,0 = 21 
y-if_(r)e< m+1 ^ J 

with integral number m. The function ip m (r,<f)) is an eigenfunction of the conserved total 
angular momentum J 3 — L 3 + S 3 — —id/d<f> + cr 3 /2 with eigenvalue j = m + 1/2. It should 
be reminded that m is not a good quantum number. This is evident from the fact that the 



7 



two components of ip m depend on the integer m in an asymmetric way. Only the eigenvalues 
j of the conserved total angular momentum J3 are physically meaningful. 
From the identities 

(T p = i(o-.r) ^-d r + ^(a-L)^j , (22) 
tr A = z(<r ■ r)a 3 A^ (23) 

one gets 

o- • (p - A)) = i(a ■ r) ( -0 r + i U3J3 - - <*A^j . (24) 



Upon using the relation 




n p im4> 

-J _ | = | ., +1U |, (25) 



one can reduce (|T8l) to 



This shows that /+ and /_ forms a one- dimensional SUSY pairs. If we take the prepotential 
W such that 

W' = 9--A4, g = m + l, (28) 
r 2 

and A ± = ±d/dr — W (as usual, the prime here denotes derivatives with respect to the 
basic variable, which is r in this case), then eqs.()26l) and (J27l) become 

A-A + f + = (E 2 -M 2 )f + , (29) 
A + A-f^ = (£ 2 -M 2 )/_, (30) 

Explicitly, the above equations read 

(~ + W' 2 ± uA f ± = {E 2 - M 2 ) f ± . (31) 

The ground state, with E 2 = M 2 , is given by one of the following two sets of equations: 

A+f?\r) = 0, f m (r) = 0, (32) 
A-f(°\r) = 0, /| 0) (r) = 0, (33) 



depending on which solution is normalizable. The solutions are generally given by 



/i 0) oc r ±(m+ ^) exp (t f drA^j . 



(34) 



To be specific, we consider the situation where m > and J drA^ > 0, so that /. 



(0) 



is 



normalizable, and = 0. The other situation can be discussed similarly. 

Eq. (|28|) relates A^ir) and W(r). This gives a way to obtain A^ that defines exactly 
solvable model. Particularly, from the Table (4.1) in one concludes that there are 
three forms of A^ giving exact solutions of the problem based on the conventional classical 
orthogonal polynomials: 

i) oscillator- like : A^{r) oc r ; 

ii) Coulomb potential-like : A^ir) oc constant ; 

iii) zero field-like : A^(r) oc 1/r . 

Case (i) corresponds simply to the well-known Landau level problem. 

Now with the discovery of the exceptional Laguerre polynomials, one can find an infinite 
family of vector potentials A^ that give the oscillator-like spectra. These are the deformed 
Landau systems. There are, however, two different types of deformed Landau systems, 
depending on whether we choose A~A + in ( J29l) to correspond to H^ + \g) in ([6]), or to 
uf\g) in (JED, or equivalently, (fl4l). 

If we choose A~ A + to correspond to H^ + \g), then from (125]) and (JTj) we get the required 
vector potential (we add a superscript to indicate its family £) 

i 



A 



m + 



uor 



2 - - w' h 



(35) 



As before, the ^ is given by ([8]) for the LI and the L2 case. These deformed systems are 
iso-spectral to the relativistic Landau system in case (i) mentioned before: the eigen-energies 
being 

StM = E ln ~M 2 = 4nu, (36) 

which is independent of g and £. The eigenfunction /+ is given by (pe, n (r; g) in ([9]): /+ oc 
<t>i,n{T; g), and /_ is obtained from /+ by ( |26l) : /_ oc A + f + . Thus both the upper and lower 
components of the eigenfunctions are related to the exceptional orthogonal polynomials in 
this case. 
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In the other choice of the vector potential, i.e., with A A + corresponding to 1%^ (<?), 
only the lower component /_ involves the new polynomials. This time (r) is determined 
by and (dU>: 

m \ ~«>r - <-=± - ®p4 : LI 

\ ujr + ^+i- : L2. 1 ^ 

(+)( 



The eigen-energies are given by £) n '(g) in f fl6|) and ([T 

m, . i o I 4 (n + g + 2£- ±) w : LI 
4t ) (^) = ^ 2 n-M 2 = , ' ( 38 ) 

[4(n + ^ + |)a; : L2. 

The eigenfunction /+ is given by (r; in ( !IB|) and (ITT]): /+ oc gr), and /_ is again 

obtained from / + by (|2T|) : /_ oc (pf^ = A^f+. Hence, /_ is related to the new exceptional 
orthogonal polynomial, as discussed in the statements below (TTTjl . Note here that in this 
case the ground state energy is not zero, and hence /_ need not vanish. Such situation is 
called broken supersymmetry in SUSY quantum mechanics 4]. 

The above discussion can be straightforwardly extended to the case of two-dimensional 
Pauli equation. This is so as the square of the Dirac equation is related to the Pauli equation. 
Specifically, we have = Hp + M 2 , where the Pauli Hamiltonian Hp is 

H P = (p- A) 2 - <r 3 (V x A) 3 . (39) 



IV. DIRAC EQUATIONS WITH NON-MINIMAL COUPLING 

The example discussed in the previous section illustrates how physical Dirac systems 
whose eigenf unctions are related to the exceptional orthogonal polynomials can be con- 
structed. The construction relies mainly on the fact that the two radial components of the 
wave function form a SUSY pair as in (1261) and (1271) . Thus, as long as a Dirac equation 
can be reduced to such a form, this construction applies and one can obtain new Dirac 
systems that involve the new polynomials. In this section we indicate three more Dirac- 
type equations which allow such construction. Unlike the the model discussed in Sect. Ill, 
these three systems involve non-minimal couplings between the fermion and the external 
electromagnetic fields. 
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A. Dirac-Pauli equation with electric field in spherical coordinates 



Consider a neutral fermion interacting with electromagnetic fields through its anomalous 
magnetic moment. The relativistic wave equation that describe such interaction is called 



the Dirac-Pauli equation [lj], ]l5| . As shown in [l5| , this equation includes as a special case 
the so-called Dirac oscillator, which has attracted much attention in recent years 16 ] . 

In this section we shall consider the case where the external field is purely electrical. 
The Dirac-Pauli equation that describes the motion of a neutral fermion of spin- 1/2 with 
mass M and an anomalous magnetic moment fi in an external electric field E is given by 
Hppty = Ety, with the Hamiltonian 

H DP = a p + ifiPa ■ E + /3M. (40) 

Here cr and (5 are the Dirac matrices which in the standard representation are given by 

where cr are the Pauli matrices. We also define \I/ = (x> where t denotes transpose, and 
both x an d v 9 are two-component spinors. Then the Dirac-Pauli equation becomes 

cr ■ (p - i//E)x = {E + M)<p, 

cr ■ (p + z/iEV = (E - M) X . (42) 

This exhibits the intrinsic SUSY structure of the system. As in the previous case, this allows 
us to construct exactly solvable Dirac-Pauli equations which are related to the exceptional 
polynomials. This we shall do below in the spherical and cylindrical coordinates. 

First let us consider central electric field E = E r (r)r. In this case, one can choose a 
complete set of observables to be {H, J 2 , J 2 ,S 2 = 3/4, A"}. Here J is the total angular 
momentum J = L + S, where L is the orbital angular momentum, and S = |S is the spin 
operator. The operator K is defined as K = 7°(S • L + 1), which commutes with both H 
and J. Explicitly, we have 

K = diag —k^j , 

k = o-L + 1. (43) 
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The common eigenstates can be written as 15 1 



here yj mj (9,<p) are the spin harmonics satisfying 



J'3?m, = ./(./ • 1)3$,,. J = -,-,..., (45) 

= '",y%„,- ki < j, (46) 

ky* mj = -ky* mj , k = ±{ ] + l -\ (47) 



and 



Using the identity (1221) . one gets 



cr • (p ± i/xE) = i(<r ■ r) ( -<9 r + ^— - ± y,E r J . (49) 
Eq.f l42|) then reduces to 

± + !L + rE^f + = (E + M)f_, (50) 

A + t + f^f. = (E-M)f + . (51) 

These two equations have the same SUSY structure as ( |26i) and ( 1271) . and so one can proceed 
as in the last section. Again, we consider the situation where k < and J dr[iE r > 0, so 
that for unbroken supersymmetry, the ground state has upper component normalizable, 
and lower component /_ = 0. The other situation can be discussed similarly. In this case, 
eq. (!28|) becomes 

W'=9--iiE t g = \k\. (52) 

Thus all the discussions following ( 128|) can be carried over with the change A$(r) — > /j,E r (r). 

Particularly, the case with fiE r oc r is just the Dirac oscillator [16|. Thus the new exactly 
solvable models associated with the Xe Laguerre polynomials constructed according to the 
procedure in Sect. Ill are the new deformed Dirac oscillators. 
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B. Dirac-Pauli equation with electric fields in cylindrical coordinates 



We now turn to 2 + 1 dimensions where the electric field is cylindrically symmetric. This 
discussion can be readily extended to the case in 3 + 1 dimensions where the electric field is 
constant along the z direction. 

The Dirac-Pauli Hamiltonian takes the form 

H DP = a ■ p + ifia 3 cr ■ E + a 3 M. (53) 

All vectors lie in the x—y plane. For E = E r (r)r, the second term in (1531) is —iix{cr-Y)<j 3 E r {r) 
(note that a 3 • <x = — a ■ 03 in this case). Comparing this term with —a ■ A in ffl~8j) and fl23l) . 
one sees that the Hamiltonians (|53|) and (|T8|) are equivalent if we make the correspondence 
fiE r (r) -H- A^(r). Therefore, the results in Sect. Ill can be directly carried over to the 
present case. 

We note here that the correspondence fiE r (r) qA^jr) (we restore the charge q) also 



underlies the duality between the Aharonov-Casher effect [17|, which is described by the 
Dirac-Pauli equation, and the well-known Aharonov-Bohm effect is| . described by the Dirac 
equation with minimal coupling. The Aharonov-Casher effect concerns the topological phase 
experienced by a neutral fermion with a magnetic moment when diffracted around a line 
of electric charge. It is an electrodynamic and quantum-mechanical dual of the Aharonov- 
Bohm effect, which gives a phase shift for a charged particle diffracting around a tube of 
magnetic flux. 

C. Two-dimensional Dirac equation with Lorentz scalar potential 

Let us consider a (1 + l)-dimensional Dirac Hamiltonian of the kind 

H = ap + f3(M + V s (x)) (54) 

where M is the mass of the fermion, p = —id/dx, a and /3 are the Dirac matrices, and V s is 
;he Lorentz scalar potential. Such model is of interest in the theory of nuclear shell model 
I9I ]. and as a model of the self-compatible field of a quark system 20]. 



This system is supersymmetric [211], as can be easily shown as follows. We represent the 
Dirac matrices by 

a = a"2, /3 — 0i. (55) 
13 



Then the Dirac equation Hip = Eip for the two-component wave function 
takes the form 

(-^-W'(x)\iP + = Eip., 

(~-W'{x)U. = E1>+. (57) 

Here W'{x) = -(V s (x) + M) and £ = E 2 . Eq. (ETJ) is now in the SUSY form, with W{x) 
playing the role of the prepotential. As such, following the procedure in Sect. Ill, we can 
obtain new exactly solvable systems related to the exceptional Laguerre polynomials. 

In fact, in this case the domain need not be confined to the half-line. Thus we can 
construct new solvable Dirac systems whose eigenfunctions are related to the exceptional 
Jacobi polynomials by simply linking W with the prepotential corresponding to exceptional 
Jacobi polynomials. 



V. FOKKER-PLANCK EQUATIONS 



Finally, we discuss briefly how the exceptional orthogonal polynomials can appear in the 
Fokker-Planck (FP) equations. In one dimension, the FP equation of the probability density 
V(x,t) is [22 1 

d 



Ot 



V(x,t) = CP(x,t), 



C 



(58) 



The functions D^'(x) and D^ 2 \x) in the FP operator C are, respectively, the drift and the 
diffusion coefficient (we consider only time-independent case). The drift coefficient represents 
the external force acting on the particle, while the diffusion coefficient accounts for the effect 
of fluctuation. Without loss of generality, in what follows we shall take = 1. The drift 
coefficient can be defined by a prepotential W(x) as D^(x) = 2W '(x) . 



The FP equation is closely related to the Schrodinger equation 



V(x,t) = e- xt e w(x U(x). 



22 



231 ] . Substituting 



(59) 
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into the FP equation, we find that <fi satisfies the Schrodinger-like equation: Hcf) = A0, where 

H = -e~ w Ce w 
d 2 



+ W'(x) 2 + W"(x). 



dx 2 

Thus satisfies the time-independent Schrodinger equation with Hamiltonian H and eigen- 
value A, and O = exp(W) is the zero mode of H: H<f) = 0. 

It is now clear that FP equations transformable to exactly solvable Schrodinger equations 
can be exactly solved. If all the eigenf unctions <p n (n = 0, 1,2, . . .) of H with eigenvalues 
A„ are solved, then the eigenf unctions V n (x) of C corresponding to the eigenvalue — A n is 
V n {x) = <po(x)(j) n (x). The stationary distribution is Vq = 0q = exp(2W) (with J Vq(x) dx = 
1), which is obviously non-negative, and is the zero mode of C: CVq = 0. Any positive 
definite initial probability density V(x, 0) can be expanded as V(x, 0) = 4>q{x) ^2 n c n 4> n (x), 
with constant coefficients c„. (n = 0, 1, . . .) 

</>nO) (0o 1 ( X )'P( X > °)) dx - ( 60 ) 
Then at any later time t, the solution of the FP equation is V(x,t) = 

4>o{ X ) J2n C n4>n{x) exp(-\ n t). 

Thus all the shape-invariant potentials in supersymmetric quantum mechanics give the 
corresponding exactly solvable FP systems. One needs only to link the prepotential W 
in the Schrodinger system with the drift potential corresponding to the drift coefficient 
£)W = 2W' in the FP system. For example, the shifted oscillator potential in quantum 
mechanics corresponds to the FP equation for the well-known Ornstein-Uhlenbeck process 



22] . 



Of interest to us here is the FP equation that corresponds to the radial oscillator po- 



tential. This FP equation describes the so-called Rayleigh process [24|. When W(x) is 
replaced by the the prepotential ©, we obtain an exactly solvable FP equation, describing 
a deformed Rayleigh process, whose eigenfunctions are given by the exceptional Laguerre 
polynomials. 

Again, as with the Dirac equation with Lorentz scalar potential, in this case one can get 
exactly solvable FP equations whose eigenfunctions are related to the exceptional Jacobi 
polynomials by simply linking W with the prepotential corresponding to exceptional Jacobi 
polynomials. 
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VI. SUMMARY 



The discovery of the exceptional Xi Laguerre and Jacobi polynomials has opened up new 
avenues in the area of mathematical physics and in classical analysis. Unlike the well-known 
classical orthogonal polynomials which start with constant terms, these new polynomials 
have lowest degree t = 1,2,..., and yet they form complete set with respect to some positive- 
definite measure. Some properties of these new polynomials have been studied, and many 
more have yet to be investigated. 

In this paper we have presented some physical models in which these new polynomials 
could play a role. For clarity of presentation, we concentrate mainly on the exceptional 
Laguerre polynomials. We show how some Dirac equations coupled minimally and non- 
minimally with external fields, and the Fokker-Planck equations can be exactly solvable with 
the exceptional Laguerre polynomials as the main part of the eigenfunctions. The systems 
presented here have enlarged the number of exactly solvable physical systems known so far. 
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